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PREFACE 


The  numerical,  steady-state,  wind-driven  circulation  study  of 
the  effect  of  a proposed  offshore  airport  on  steady-state  circulation 
in  Lake  Erie  near  Cleveland,  Ohio,  was  sponsored  by  the  Lake  Erie  Regional 
Transportation  Authority  (LERTA),  Cleveland,  Ohio,  as  a part  of  the 
model  feasibility  investigation  being  conducted  at  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES).  The  WES  investigation,  Task  17  of 
the  LERTA  investigation,  is  a portion  of  the  second-phase  airport  feasi- 
bility study  undertaken  by  LERTA  to  evaluate  proposed  airport  sites,  one 
of  which  is  in  Lake  Erie  near  Cleveland. 

This  report  and  the  numerical  analyses  were  prepared  and  conducted 
by  Dr.  Donald  L.  Durham  and  Mr.  H.  Lee  Butler  of  the  Wave  Dynamics 
Division  (WDD),  WES,  under  the  general  supervision  of  Dr.  R.  W.  Whalin, 
Chief,  WDD,  and  Mr.  H.  B.  Simmons,  Chief,  Hydraulics  Laboratory. 

Assisting  in  this  study  were  Messrs.  D.  G.  Outlaw  and  C.  W.  Coe,  WDD, 
and  Dr.  Donald  C.  Raney,  who  was  working  with  WDD  while  on  loan  from 
AME  Department,  University  of  Alabama  through  the  Intergovernmental 
Personnel  Exchange  Program. 

The  Director  of  WES  during  the  conduct  of  this  investigation  and 
the  preparation  and  publication  of  this  report  was  COL  G.  H.  Hilt,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 

Multiply By To  Obtain 


inches 

2.54 

centimeters 

inches  per  second 

2.54 

centimeters  per  second 

inches2  per  second 

6.4516 

centimeters  per  second 

pounds 

4. 448xl05 

dynes 

feet  per  second 

0.3048 

meters  per  second 

miles 

1.6093 

kilometers  per  hour 

miles  per  hour 

1.6093 

kilometers  per  hour 

PART  I : INTRODUCTION 


1.  The  Lake  Erie  Regional  Transportational  Authority  (LERTA) 
is  conducting  a feasibility  and  site  selection  study  for  a major  hub 
airport  in  the  Cleveland  Service  area.  One  of  the  possible  sites  being 
evaluated  is  an  offshore  site  in  Lake  Erie  near  Cleveland,  Ohio.  As 
a part  of  the  feasibility  analysis  of  an  offshore  site,  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  conducted  a model  feasibility 
investigation  and  is  performing  numerical  model  studies.  The  results 
of  these  efforts  are  being  published  in  a series  of  reports  under  the 
general  title  "Lake  Erie  International  Jetport  Model  Feasibility  Investi- 
gation." This  miscellaneous  paper  presents  preliminary  results  of  the 
numerical  steady-state,  wind-driven  circulation  analyses  for  a lake  of 
constant  density.  These  data  are  the  results  of  a parametric  study 
(wind-direction  and  speed)  of  the  steady-state,  wind-driven  circulation 
with  and  without  a jetport  island  in  Lake  Erie  near  Cleveland,  Ohio. 


PART  II:  HYDRODYNAMIC  MODEL 


2.  From  earlier  efforts  of  WES  in  the  model  feasibility  study, 
the  numerical  model  which  is  used  in  this  study  was  selected1  from 

2 

several  numerical  steady-state,  wind-driven  circulation  models  and  applied 
to  test  cases  of  a jetport  island  offshore  of  Cleveland  in  Lake  Erie. 

i 

Details  of  the  hydrodynamic  model  and  the  development  of  the  finite 
difference  scheme  for  this  model  are  presented  in  References  1 and  3. 

A cursory  statement  of  the  hydrodynamic  model  and  assumptions  are  presented 
here  as  a brief  summary  of  the  model. 

3.  The  governing  partial  differential  equation  for  the  integrated 
stream  function  i//*  is  the  following  Poisson  equation: 

V2^  + o2<J>y  = 03  Cl) 

where  o^,  o2,  and  o 3 are  functions  of  local  depth  and  bottom  slope  of 
Lake  Erie.  Also,  a3  is  a function  of  wind  stress.  Detailed  definitions 
of  these  functions  and  the  complete  development  of  the  equations  can  be 
found  in  References  1,  2,  and  3.  The  integrated  stream  function  ij; 
is  defined  by  U = 3T/3y  = i )jv  and  V = - 3i|//3x  = where  U is  the 

vertically  integrated  velocity  component  in  the  direction  of  the  x axis 
and  V is  the  vertically  integrated  velocity  component  in  the  direction 
of  the  y axis.  The  boundary  condition  for  this  problem  is  that  \p  be 
constant  along  the  shore. 

4.  To  numerically  solve  the  stream  function  problem,  the  governing  , 

partial  differential  equation  (Equation  1)  is  approximated  by  finite 
differences^  at  a discrete  number  of  grid  points  within  the  flow  field 

i 

of  Lake  Erie.  Plates  1 and  2 depict  the  discretized  grids  for  the 

' 

entire  Lake  (two-mile  grid  spacing)  and  for  the  12-mile  by  12-mile 
nearshore  area  (0.25  mile  grid  spacing)  near  Cleveland,  Ohio.  The  finite 
difference  approximation  of  Equation  1 at  a finite  number  of  grid  points 
produces  a finite  set  of  simultaneous  equations  governing  . This 
set  of  simultaneous  equations  is  solved  iteratively  by  a successive 

* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 
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over  relaxation  method  to  determine  at  each  grid  point  a value  of  \p  , • 

which  satisfies  the  partial  differential  equation  (Eq  1)  and  the  specified 
boundary  condition.  Knowing  the  stream  function,  the  pressure,  surface 
displacement,  and  vertically  integrated  velocity  components  can  be 
obtained. 

5.  For  the  u and  v components  of  horizontal  velocity,  Gedney 
developed3  analytical  expressions,  which  are  functions  of  lake  depth, 
wind  stress,  gradient  of  surface  displacement,  and  vertical  distance, 

z , from  the  still  water  level.  Thus,  by  knowing  < 1/  the  horizontal 
velocity  can  be  determined  at  any  depth  in  the  water  column.  Using 
Gedney' s expressions  of  u and  v and  integrating  in  closed  form  the 
continuity  equation,  an  analytical  expression  for  the  vertical  component, 
w , of  velocity  can  be  obtained  as  a function  of  wind  stress,  lake 
depth,  surface  displacement,  and  z . Therefore,  the  complete  three- 

! 

dimensional  flow  field  can  be  determined. 

6.  For  this  numerical  study  of  the  effects  of  a proposed  jetport 
island  on  the  steady-state,  wind-driven  circulation  for  well-mixed  conditions 
in  Lake  Erie,  the  bottom  topography,  wind  stress,  and  various  other  character- 
istic parameters  must  be  specified  as  input  to  the  study.  Discrete  values 

of  the  bottom  topography  of  Lake  Erie  were  obtained  for  each  grid  point  of 
the  two  grids  shown  on  Plates  1 and  2.  For  the  entire  Lake  (Plate  1),  discrete 
values  of  depths  were  obtained  from  Lake  Survey  Center  Charts  as  published 
by  the  U.  S.  Army  Corps  of  Engineers.  Reference  3 outlines  the  procedure 
with  which  these  values  were  read  from  charts,  interpolated  and  smoothed. 

Contours  of  the  resultant,  discretized  bottom  topography  for  the  full 
Lake  are  depicted  in  Plate  3.  For  the  finer  nearshore  grid  (Plate  2), 
discrete  values  for  depth  were  obtained  by  interpolating  over  the  larger 
full  lake  grid.  Near  the  shore  (within  5 miles),  additional  information 
* I-  as  to  the  bottom  topography  and  any  abrupt  changes  in  depth  was  obtained  from 

other  available  maps.  Contours  of  the  resultant,  discretized  bottom  topography 
for  the  nearshore  gird  are  depicted  in  Plate  4. 

7.  The  wind  data  used  in  this  s.tudy  of  the  steady-state,  wind- 
driven  circulation  arc  based  on  an  analysis  (Reference  5)  of  wind 


i 


conditions  offshore  of  Cleveland,  Ohio.  Ten  years  of  recorded  wind  data 
from  Burke  Lakefront  Airport  and  Cleveland  Hopkins  Airport  were  analyzed 
by  A.  H.  Glen  § Associates5  to  obtain  average  monthly  and  annual  wind  speed 
direction  regimes  and  percentage  frequency  of  occurrence  for  two  off- 
shore stations  near  the  proposed  jetport  site.  For  this  steady-state 
study,  the  lake  water  is  assumed  to  be  isothermal,  which  implies  constant 
density.  This  condition  in  Lake  Erie  generally  exists5  from  late 
September  through  May.  On  the  basis  of  the  wind  study,5  two  wind  magnitudes 
(17  mph  and  35  mph)  were  chosen  as  representative  of  average  and  extreme 
wind  conditions  during  the  time  period  from  late  September  through  May. 
Winter  winds  in  the  Cleveland  area  most  frequently  are  from  south  to 
west-southwest.  However,  the  wind  direction  for  this  study  was  considered 
a parameter  and  for  each  of  the  two  wind  magnitudes  was  varied  in  22.5° 
increments  from  0°  to  360°.  Therefore,  16  wind  directions  for  each  wind 
magnitude  were  considered  in  the  numerical  study.  For  each  wind  field 

/v 

(magnitude  and  direction),  the  wind  shear  stress,  t,  is  determined  by 

. , 2 

| T | = 0.00273  paWa  (2) 


where 

p„  = density  of  air, 

a 

W = wind  velocity  in  m/sec, 

and  0.00273  is  the  drag  coefficient  by  Wilson's0  formulation. 

8.  The  inflow  and  outflow  of  the  major  rivers  for  Lake  Erie  are 
considered  in  specifying  the  boundary  value  of  the  stream  fuction 
along  the  shore.  The  major  rivers  which  account  for  the  through-flow 
of  Lake  Erie  are  the  Detroit,  Maumee,  Sandusky,  and  Niagara.  Although 
the  inflow  of  the  Cuyahoga  River  has  only  a small  local  effect  on  the 
Lake  circulation,  this  effect  is  of  interest  to  this  study.  Thus,  the 
Cuyahoga  River  inflow  is  also  considered  in  specifying  the  boundary 
condition.  The  boundary  condition  is  that  the  stream  function 

along  the  shore  is  a specified  function,  ip  = ip(s) , of  the  shoreline 
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To  account  for  river  discharge  in  specifying  the  stream  function 

along  some  portion  of  the  shore  is  arbitrarily  set  to  a constant 
(e.g.  0.0)  with  changes  in  the  specified  value  of  the  stream 

function  occurring  at  the  mouth  of  major  rivers  and  representing 
the  volume  flow  rate  of  river  inflow  or  outflow.  Rivers  and  average 
volume  flow  rates  used  in  this  study  are  present  in  Table  1. 

2 

9.  In  the  non-dimensional  formulation  of  the  equations  of  motion 
various  characteristic  parameters  are  defined  and  are  used  in  scaling  the 
non-dimensional  results  of  the  numerical  study.  A list  of  the  pertinent 
parameters  and  the  values  used  in  this  study  is  presented  in  Table  2.  In 
addition  to  these  characteristic  parameters,  the  frictional  depth  D 
vertical  eddy  viscosity  Ay  , and  the  magnitude  |t|  of  the  wind  shear 
stress  are  dependent  on  the  magnitude  of  the  forcing  wind.  Calculated 

7 ~ 1 2 i 

values  of  75  ft  and  115  ft  for  D,  27  cnT  sec  and  60  cn  sec'1  for  Ay, 
and  1.9  dyne  cm’2  and  S.55  dyre  cm’2  for  J t | were  used  respectively  for 
wind  speeds  of  17  mph  ar.d  35  mpn. 


! 

I 

i 
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PART  III:  RESULTS  OF  NUMERICAL  STUDY 


10.  The  results  of  the  numerical  steady-state,  wind-driven  circulation 
analysis  of  the  effects  of  a proposed  jetport  island  on  circulation  in 
Lake  Erie  are  presented  graphically  as  Appendices  B and  C.  The  steady- 
state  wind  conditions  considered  in  the  analysis  are  16  wind  directions 
at  22.5°  increments  from  0 to  360°  with  wind  speeds  of  17  mph  and  35  mph. 

For  9 of  the  16  wind  directions  (direction  from  which  wind  is  blowing), 
the  results  have  been  consolidated  on  data  plates  and  are  presented  as 
Appendix  B.  These  results  in  Appendix  B are  grouped  by  wind  speed 
(17  mph  results  first)  and  direction  in  the  following  order: 


East 

North 

Northwest 

West 

West  southwest 
Southwest 
South  southwest 
South 

South  southeast 


For  each  of  these  directions  and  magnitudes,  ten  data  plates  are  presented. 
The  wind  speed  and  direction  are  given  on  each  data  plate,  e.g.  17  mph  E 
(17  miles  per  hour  from  the  east).  The  data  on  each  plate  for  a specific 
speed  and  direction  are  as  follow: 

Plate  1:  Stream  function  contours  and  horizontal  velocities  for 

Lake  Erie. 

Note:  The  stream  function  contours  are  vertically 

integrated.  The  horizontal  velocities  are  at  the 
lake  surface  and  at  the  5-ft  depth. 

Plate  2:  Horizontal  velocities  at  the  10- , 20- , 40- , and  60-ft 

depths  (for  the  full  lake). 

Plate  3:  Vertical  velocities  at  10-,  20-,  40-,  and  60-ft  depths 

(for  the  full  lake). 

Note:  The  stream  function  contours  are  vertically  inte- 

grated. The  full  lake  solution  has  been  used  as  an 
offshore  boundary  condition  for  the  nearshore  grid  solution. 


9 


1 


Plate  5:  Nearshore  horizontal  velocities  at  0-,  10- , and  15-ft 

depths. 

Plate  6:  Nearshore  horizontal  velocities  at  20-,  30- , and  40-ft 

depths.  ' 

Plate  7:  Nearshore  vertical  velocities  at  the  10-,  20- , and  30-ft 

depths. 

Note:  The  length  scale  for  the  velocity  magnitude  has 

been  increased  by  the  amplification  factor  shown. 

Plate  8:  Nearshore  contours  of  differences  in  velocity  components 

or  magnitude  at  a depth  of  0 ft. 

Note:  The  direction  of  the  U component  of  the 

velocity  is  along  the  grid  line  approximately  parrel  lei 
to  the  shoreline  shown  on  the  plot.  The  direction  of  the 


V component  of  the  velocity  is  along  the  grid  line 
approximately  perpendicular  to  the  shoreline.  The  velo- 
city difference  was  obtained  by  subtracting  the  velocity 
calculated  with  a jetport  island  in  the  Lake  from  the 
velocity  calculated  without  a jetport.  The  velocity 
magnitude,  M , difference  was  obtained  in  the  same 
manner  but  without  regard  to  direction. 


Plate  9: " Nearshore  contours  of  difference  in  velocity  components 
or  magnitude  at  a depth  of  15  ft. 

Note:  Same  comments  as  for  Plate  8. 

Plate  10:  Nearshore  contours  of  differences  in  velocity  components 

or  magnitude  at  a depth  of  30  ft. 

Note:  Same  comments  as  for  Plate  8. 

11.  For  the  remaining  seven  directions,  consolidated  data  plates 
are  not  prepared;  however,  the  same  information  is  available  from  plots 
of  the  original  data.  These  plots,  reduced  in  size,  are  presented  as 
Appendix  C.  In  a similar  manner  as  Appendix  B,  the  data  in  Appendix  C 
are  grouped  (1)  with  the  17  mph  wind  speed  plots  for  the  east 
northeast,  northeast,  north  northeast,  north  northwest,  west  northeast, 
southeast,  and  east  southeast  directions  and  (2)  with  plots  of  the 
35  mph  wind  speed  for  the  same  directions.  Forty  plots  are  included 
for  each  of  the  wind  speeds  and  directions.  Table  3 cross  references 
plots  in  Appendix  C with  data  plates  and  associated  comments  of  Appendix  B. 
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12.  The  numerical  model^,  which  was  developed  at  Case  Western 
University  to  compute  steady-state,  wind-driven  currents  for  the  full 
lake  and  nearshore  region  of  Lake  Erie,  consisted  of  eleven  separate 
codes.  WES  combined  these  eleven  codes  into  two  jobstreams  for  efficient 
execution  on  a GE  635  computer.  In  addition  to  converting  these  codes, 
WES  added  contouring  and  plotting  routines.  The  first  jobstream  obtains 
a solution  for  the  full  lake.  The  second  jobstream  obtains  a solution 
for  the  nearshore  region  near  Cleveland  Harbor.  A description  of  each 
activity,  input/output  requirements  and  definitions,  flowchart,  and  deck 
setup  are  provided  for  each  jobstream  in  Appendix  D.  Appendix  D is 
bound  and  published  separately. 


PART  IV:  CONCLUSIONS 


13.  Based  on  the  numerical  results  of  the  parametric  study 
(wind  direction  and  speed)  of  the  effects  of  a jetport  island  on 
the  steady-state,  wind-driven  circulation  for  constant  density  condi- 
tions in  Lake  Erie  near  Cleveland  Ohio,  it  is  concluded  that: 

a.  For  all  wind  directions  and  magnitudes,  the  effect  of  the 
jetport  island  on  lake  cir-jlation  is  not  appreciable 
except  within  2 miles  of  the  island. 

b.  Perturbations  in  the  velocity  field  for  the  17  mph  and 
35  mph  wind  fields  are  very  similar  with  two  exceptions 
for  the  35  mph  wind  fields.  These  exceptions  are 

(1)  a slight  increase  in  surface  area  over  which  the 
velocity  perturbations  are  found  and  (2)  an  increase 
in  intensity  of  these  perturbations  very  close  to  the 
island. 

c.  Two  windows  of  wind  directions  (WNW  through  WSW  and  ESE 
through  ENE)  are  characterized  by  the  largest  perturbations 
(both  magnitude  and  areal  extent)  in  the  velocity  field 
due  to  the  effects  of  the  jetport  island. 

d.  Wind  from  the  N or  S directions  generate  circulation 
patterns  which  possess  minimum  effects  associated  with 
the  jetport  island. 

e.  For  the  above  two  wind-direction  windows,  the  maximum  effect 
(relative  to  other  wind-directions)  of  the  island  on  the 
steady-state  circulation  patterns  in  the  nearshore  region 

is  due  to  the  island  location  being  in  a region  of  large 
net  mass  flux.  Thus,  the  island's  presence  in  such  a region 
produces  substantial  blockage  of  the  net  flow.  Comparisons 
of  plots  of  streamlines  for  each  wind  condition  vividly 
demonstrates  this  result  (e.g.  streamlines  for  17  mph  W 
results  versus  streamlines  for  17  mph  N results). 

f.  During  the  winter  winds  or  well-mixed  lake  conditions 
(September  through  May),  winds  occur  only  8 percent 
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of  the  time  from  the  ENE-ESE  window  and  22  percent  of  the 
time  from  the  KSW-WNW  window5.  Annual  frequency  of  occurrences 
for  these  windows  are  8 and  19  percent  respectively.'  The  most 
frequently  occurring  wind  directions  are  contained  in  a window 
from  SSE  to  SW  with  an  annual  frequency  of  occurrence  of 
42  percent  and  42  percent  for  winter  wind  conditions.  Thus, 
for  the  well-mixed  lake  conditions,  the  two  wind-direction 
windows  which  are  characterized  by  maximum  island  effects  do 
not  possess  the  most  frequently  occurring  wind  directions  and 
have  low  frequencies  of  occurrence.  Therefore,  the  maximum 
effects  of  the  island  on  the  steady-state  circulation  patterns 
will  occur  much  less  frequently  than  the  smaller  effects 
associated  with  other  wind  directions. 

For  all  wind  directions,  maximum  horizontal  velocity  values 
at  the  lake  surface  between  the  jetport  and  the  shoreline  are 
1.5  ft/sec  for  the  17  mph  wind  speed  and  2.5  ft/sec  for  the 
35  mph  wind  speed. 

The  jetport  effects  in  the  horizontal  velocity  region  near 
the  shoreline  in  the  vicinity  of  Cleveland,  Ohio  are  estimated 
to  be  less  than  0.1  ft/sec  for  all  studied  wind  fields. 

Upwelling  and  downwelling  along  the  shoreline  in  the  vicinity 
of  Cleveland,  Ohio  are  affected  very  little,  if  any,  with  no 
changes  in  location  along  the  shoreline  of  the  areas  of 
upwelling  and  downwelling  and  only  slight  changes  for  a few 
wind  directions  (i.e.,W)  in  magnitude  of  vertical  velocities 
in  these  areas. 

Around  the  perimeter  of  the  island,  areas  of  upwelling 
and  downwelling  are  evident  with  the  location  of  these 
areas  depending  on  wind  direction.  The  magnitude  of 
vertical  velocities  are  extremely  small  and  of  the  order 
of  10"3  to  10'4  ft/sec. 
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TABLE  1 


Average  River  Discharges  (cfs) 


Cuyahoga 

Detroit 

Maumee 

Sandusky 

Niagara 


850 

160,000 
4,825 
1,045  #1 
-166,720 


TABLE  2 

Characteristic  Parameters 


H Characteristic  depth 

m 

L Characteristic  length 

Trefl  Reference  wind  stress 

f Coriolis  parameter 

u _ Reference  velocity 


18.2  m (60  ft) 

128.75  km  (80  miles) 
-2 

1 dyne  cm 
10'4  sec'1 
5.47  cm  sec'1 
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TABLE  3 


Plot  Cross-Reference 


Data  Plate  No. 


Plot  No. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


1-3 

4-7 

8-11 

12,13 

14-19 

20-25 

26-31 

32-34* 

35-37* 

38-40* 


i 

i 

! 


j 

* The  velocity  component  (U,V)  or  magnitude  (M),  depth  (p),  wind  speed, 
and  direction  arc  labeled  on  plots  32-40  as  U,  D=0,  W= 1 7I-NH  to  indicate 
that  the  plot  is  for  the  U velocity  component;  the  depth  is  0-ft;  the 
wind  speed  is  17  mph;  and  the  direction  is  cast  northeast,  for  example. 
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SYMBOLS  CONTOURS.  FT 

1 0.000 

2 20.000 

3 30 .000 

4 40.000 

6 46.000 

• 60.000 

T 65.000 

0 00.000 
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CLEVELAND  NEARSHORE  AREA 


NOTATION 


Ay  vertical  eddy  viscosity  coefficient 
D frictional  depth  and  lake  depth 
characteristic  depth 
L characteristic  length 

u horizontal  velocity  in  x direction  at  given  depth 

U vertically  integrated  velocity  component  in  the  direction  of 

the  x axis 

U _ reference  velocity 
ref 

v horizontal  velocity  in  y direction  at  given  depth 

V vertically  integrated  velocity  component  in  the  direction  of 

" the  y axis 

w vertical  component  of  velocity 
Wa  wind  velocity  in  m/sec 
x,y,z  Cartesian  coordinates 
Aip  change  in  stream  function 
Pa  density  of  air 

Oj.a 2 function  of  local  depth  and  bottom  slope  of  Lake  Erie 
0 3 function  of  local  depth,  bottom  slope,  and  wind  stress  of  Lake  Erie 


Sip/  3x 
3ip/3y 


wind  shear  stress 
magnitude  of  wind  shear  stress 
reference  wind  stress 
integrated  stream  function 

boundary  value  of  stream  function  specified  along  the 
shore  line 

first  derivative  of  ip  with  respect  to  x 
first  derivative  of  ip  with  respect  to  y 
Laplacian  of  ip 
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